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ABSTRACT 

Sub-TeV gamma-ray emission from the North- West rim of the supernova remnant RX J0852. 0—4622 
was detected with the CANGAROO-II telescope and recently confirmed by the H.E.S.S. group. In 
addition, the H.E.S.S. data revealed a very wide (up to two degrees in diameter), shell-like profile of 
the gamma-ray emission. We have carried out CANGAROO-III observations in January and February 
2005 with three telescopes and show here the results of three-fold coincidence data. We confirm the 
H.E.S.S. results about the morphology and the energy spectrum, and find the energy spectrum in the 
NW-rim is consistent with that of the whole remnant. 

Subject headings: gamma rays: observation — supernova remnant: individual (RX J0852. 0—4622) 



1. INTRODUCTION 

RX J0852.0-4622 (G266.2-1.2) is one of the few 
supernova remnants (SNRs) from which strong non- 
therma,l X-ray emission has b een detected llAschenbachI 
ITqM iTsunemi et all l200(it ISlane et all 12001 ! A 
distance of ~200pc and a n age of ^700 years 
f\schenbach et al. 1999; Tsu nemi et al.l 120001) would 
make RX J08 52.0-4622 one of the closest young SNR 
(although see ISlane et al.l 120011 for a dissenting view) . 
The non-thermal emission suggests particle acceleration 
to ^lOOTeV, with the possibility of TeV gamma-ray 
emission from inverse Compton scattering off electrons, 
and/or 7r° decay from proton interactions, with clear im- 
plications for the origin of cosmic rays in our galaxy 
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ijVnlk et a1jl200l . 

Sub-TeV gamma-rays from the northwest rim (NW 
rim) of RX J0852.0-4622 were detected with the 10 m 
CANGAROO- II imaging atmosphe ric Cherenkov tele- 
scope (lACT) (|Katagiri et al.ll2005t) . Gamma-ray emis- 
sion was recently confi rmed by the H.E.S.S. group 
ijAharonian et aljl2005b() . who mapped the spatial distri- 
bution of gamma-ray emission over the whole remnant, 
using data from at least two of four telescopes in their 
array. 

RX J0852. 0—4622 is the seco nd spatially resolved 
SNR at TeV gam ma-ray energies (jEnomoto et al.ll2002l: 
lAharonian eTaI]l200^ . The CANGAROO-II result in- 

dicated a soft energy spectrum oc E ' however, in 
contrast the H.E.S.S. group reported a harder spectrum 
^-2.i±o.i±o.2^ rpj_^g gamma-ray fluxes also appeared dis- 
crepant: at 1 TeV the CANGAROO-II flux level was 15% 
of the Crab nebula, whereas the H.E.S.S. result was 100% 
of the Crab. The clear differences between these two re- 
ports are; 

• CANGAROO-II was only sensitive to within 0.4° 
of the X-ray maximum i n the N W rim i j A schenbachl 
1199^ ITsunemi et aTl 120001: ISla.ne et abl 120011^ 
whereas H.E.S.S. was sensitive to the whole rem- 
nant. 

• the CANGAROO-II observation used a single 
telescope whereas the H.E.S.S. observations used 
stereo data from two or more telescopes. 

• the CANGAROO-II flux was derived based on the 
subtraction of a background derived from the sur- 
rounding area which included part of the remnant, 
whereas the H.E.S.S. background was derived from 
an off-source region outside the remnant. The 
CANGAROO-II result, therefore, was dependent 
on the spatial distribution of gamma-rays inside 
the SNR. 

The H.E.S.S. result (from stereo observations) showed 
very extended emission, up to one degree in radius from 
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the center of the remnant. In this paper we present the 
results of CANGAROO-III stereoscopic observations and 
determine how much of the apparent discrepancies be- 
tween these two earher resuhs arise from these differ- 
ences. 

2. CANGAROO-III STEREOSCOPIC SYSTEM 

The use of imaging atmospheric Cherenkov telescopes 
(lACTs) was established with the statistically unassail- 
able detection o f the Crab nebula at TeV energies by the 
Whipple group (jWeekes et alJIlQSfll) . This technique en- 
ables TeV gamma-rays to be selected from the huge back- 
ground of cosmic rays with the use of the "Hillas param- 
eters" of the shower images (Hillas 1985). Stereoscopic 
observations, which allow the signal-to-noise ratio to be 
signific antly improved, were p ioneered by the HEGRA 
group ijAharonian et a n IT99I . The H.E.S.S. group 
has recently reported the detection of faint gamma-ray 
sources with an angular resolut ion as fine as a few arc- 
minutes ( Aharonia n et alJl2005D . 

CANGAROO-III is one of two major lACTs located 
in the southern hemisphere. The CANGAROO-III 
stereoscopic system consists of four imaging atmospheric 
Cherenkov telescopes located near Woomera, South Aus- 
tralia (31°S, 137°E). Each telescope has a reflector of 
10 m diameter. Each reflector consists of 114 segmented 
spherical mirrors (80 cm in diamete r with a radius of cur- 
vature of 16.4 m) made of FRP l|Kawachi et alJ I2001D 
mounted on a parabolic frame {f /d=Q.77, i.e., a focal 
length of 8m). The total light collection area is 57.3 m^. 
The first telescope, Tl, which was the CANGAROO-II 
telescope, is not presently in use due to its smaller field 
of view and higher energy threshold. The second, third, 
and fourth telescopes (T2, T3, and T4) were used for 
the observations described here. The camera systems for 
T2, T3, and T4 are identical and their details are given 
in iKabuki et al.l (j2p03). The telescopes are located at 
the east (Tl), west (T2), south (T3) and north (T4) cor- 
ners of a diamond with sides of ~100m (|Enomoto et al.l 
I2002bft . 

3. OBSERVATIONS 

The observations were carried out in the period from 
2005 January 16 to February 23 using "wobble mode" in 
which the pointing position of each telescope was shifted 
in declination between ±0.5 degree from the center of 
the remnant ( RA, dec = 133.00° , -46.37° J2000) ev- 
ery 20 minutes l)Daum et alJll997l) . Data were recorded 
for T2, T3 and T4 when more than four photomulti- 
plier (PMT) signals exceeded 7.6 photoelectrons (p .e.) in 
any t wo telescopes ("Global Trigger") (Nishiiima et al.l 
|2QQ3). The typical trigger rate was 9 Hz for three- fold 
coincidence. Each night was divided into two or three pe- 
riods, i.e., ON-OFF, OFF-ON-OFF, or OFF-ON obser- 
vations. ON-source observations were timed to contain 
the meridian passage of the target. On average the OFF 
source regions were located with an offset in RA of 4-30° 
or —30° from the center of the remnant. The OFF-source 
observations were also made in wobble mode. The total 
observation time was 1736 and 1408 min, for ON and 
OFF observations, respectively. Next we required the 
images in all three telescopes to have clusters of at least 
five adjacent pixels exceeding a 5 p.e. threshold (three- 
fold coincidence) . The event rate was reduced to ^^6 Hz 



by this criterion. Looking at the time dependence of 
these rates, we can remove data taken in cloudy condi- 
tions. This procedure is the same as the " cloud cut " used 
in the CANGAROO-II analysis (Enomotoera0|2002). 
We also rejected data taken at elevation angles less than 
60°. In total, 1081 min. data survived these cuts for ON 
and 1031 min. for OFF, with a mean elevation angle of 
70.3°. The size of this supernova remnant is large (2 de- 
gree) even compared to the field of view (FOV) of our 
observation (4 degree). One of the reason why we took 
"wobble" observation is to enlarge the effective FOV, the 
other is to average the responses of individual pixels. We, 
therefore, took LONG OFF source run of "wobble" mode 
for background subtractions in the later analysis. 

Prior to these observations, we observed this SNR 
from 2004 January to February with two telescopes 
(T2 and T3). and pr eliminary results can be found in 
iTanimori et al.l l)2005(l . The statistics of the OFF-source 
runs were insufficient, only 1/3 of the ON-source runs. 
Therefore, we describe that data set in the Appendix 
and concentrate on the description of 2005 observations 
in the body of this paper. 

The light collecting efficiencies, including the reflectiv- 
ity of the segmented mirrors, the light guides, and the 
quantum efRciencies of photomultiplicr tubes were mon- 
itored by a muon-ring analysis (Enomoto et al. 2006). 
The light yield per unit arc-length is approximately pro- 
portional to the light collecting efficiencies. The ratios 
of these at the observation period with respect to the 
mirror production times (i.e., deterioration factors) were 
estimated to be 45, 55, and 73% for T2, T3, and T4, re- 
spectively. The measurement errors are considered to be 
at less than the 5% level. These values were checked ana- 
lyzing Crab data which were obtained in 2004 November, 
which are described in the next section. The deteriora- 
tions were mostly due to dirt and dust settling on the 
mirrors. We cleaned the mirrors with water in October 
2005 and the partial improvement (a factor of 1.3-1.4) of 
the light collecting efficiencies were observed. 

4. ANALYSIS 

The analys is procedures us e d wer e identical with those 
described in lEnomoto et alJ l)2006j) for the Vela pulsar 
wind nebula, which is approximately 0.5 degree south 
of the Vela Pulsar. A s a full instrumental description 
was given in Enomoto et al. (200^, we omit a detailed 
discussion here. 

At first, the Hillas parameters l)Hillas II1985D were cal- 
culated for the three telescopes' images. The gamma- 
ray incidence directions were adjusted by minimizing the 
sum of squared widths (weighted by the photon yield) of 
the three images seen from the assumed position (fitting 
parameter). 

In order to derive the gamma-ray likeliness, we used an 
established mathematical method, the Fisher Discrimi- 
nant (hereafter FD in short) (Fisher 1936,1 . When we 
use a multi-parameter set such as: 

P = (W2, W3, WA, L2, L3, LA), 

where W2, Wi, WA, L2, L3, LA are energy corrected 
widths and lengths for the T2, T3, and T4 and assume 
that a linear combination of 

FD^d- P, 
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Fig. 1. — Fisher Discriminant (FD) distributions. The black 
points with the error bars are those for events inside a one degree 
circle centered on the SNR. The green histogram is for OFF-source 
data. The normalization of the latter was determined by the fit- 
ting procedure described in the text. The blue points show the 
background subtracted signals and the red histogram is the best 
fit for gamma-rays from Monte-Carlo simulations. 



provides the best separation between signal and back- 
ground, then the set of linear coefficients (a) should be 
uniquely determined as 

Psig — f^BG 

a = , 

Esig + Ebg 

where /I is a vector of the mean value of P for each sam- 
ple and E is their correlation matrix. This method is 
widely used in the high-energy experiments such as the 
B — factory in order to separ ate spherical events from 
jet-like ones ijAbe et all 20011) . We previousl y used it 
in Vela Pulsar analysis I Enomoto et alJ f2006j) to sepa- 
rate "sharp" (gamma-ray-like) images from "smeared" 
(background) ones. The values of jlsig, /Ibg, Esig, and 
Ebg can be calculated from the Monte-Carlo and obser- 
vational data (OFF-source runs), respectively. 

FD is a linear combination of multi-parameters. We, 
therefore, should select parameters which can be linearly 
added. The widths and lengths of T2, T3, and T4 are 
second order cumulative moments of shower images and 
thus their linear combination is a reasonable assump- 
tion. In order to keep the mean of the Hillas parameters 
(which are energy dependent) at the same values in vari- 
ous energies, we corrected the energy dependence of each 
parameter subtracting best-fitted second-order polyno- 
mials for the Monte-Carlo gamma-rays, i.e., the means 
of the corrected Hillas parameters are all set to zero for 
gamma-rays, independent of energy. This procedure en- 
sures the mean position of the FD to be located at zero 
for gamma-rays. Also a linear combination of six ran- 
dom parameters ensures the distribution has an approxi- 
mately Gaussian form (see the red histogram in Fig. ^ of 
Monte-Carlo gamma-ray events). The contributions of 
Hillas parameters to the photon-proton separation can 
be expressed by ai-Pi, which were (0.06, 0.25, 0.48, 0.28, 
0.31, 0.16) for W2, W3, WA, L2, L3, M, respectively. 



Fig. 2. — The energy resolution (AE/E) estimated by Monte- 
Carlo simulations. The red points are obtained before the "edge 
cut" and the black points after it. The curve is the best- fit function 
described in the text. 



The main differences of three telescopes are the different 
spot sizes of mirrors. As was expected, the width of T4 
gave the largest contribution which reflects the respective 
mirror qualities. The FD can be understood as an ap- 
proximate representation of "mean scaled widtK^ which 
was energy corrected and linearly corrected by the other 
Hillas parameters. The advantage of this method is that 
it can be calculated fully automatically in the analysis 
procedure, i.e., any lACT can reproduce it. 

Here, we rejected events with hits in the outermost 
layer of the cameras ("edge cut"). These rejected events 
cause finite deformations especially in the length distri- 
bution which results in deformations of the FD. This 
cut w as introduced after the analysis in lEnomoto et afl 
j200^ and was found to improve the energy resolu- 
tion significantly: Ai?/_E was estimated using Monte- 
Carlo simulations and is plotted in Fig. [3 The red 
points are those before the "edge cut" and the black 
after. A dramatic improvement in energy resolution is 
obtained. It was parametrized by the function l^E/ E = 
TJl^Ejl TeV + 2% (the solid curve in Fig.[3). Here, 
we did not apply either the "core dis tance" or "sl ant dis- 
tance" corrections described in H ofmann et al.l (|2000). 
The gamma-ray energies were simply estimated by the 
total number of photoelectrons observed. Although this 
cut significantly reduces gamma-ray acceptance, espe- 
cially in the higher energy region, we decided to use it 
in order to safely determine the gamma-ray energy spec- 
trum in this analysis. The improvement of acceptance 
may be considered in more detail elsewhere in the fu- 
ture. 

Since we have FD distributions for OFF-source data 
and the Monte-Carlo gamma-ray events, we can assume 
these are background and signal behaviors. We, there- 
fore, can fit the FD distribution of ON (the black points 
in Fig.^ with the above emulated signal and real back- 
ground functions, to derive the number of signal events 
passing the selection criteria. With this fit, we can deter- 



4 



Enomoto et al. 




Fig. 3. — The 6^ plot for the 2004 November Crab nebula data. 



mine the gamma-ray excess without any positional sub- 
tractions, i.e., we can analyze all events in the FOV and 
not only the events which comes from the direction of 
the SNR. This is a two-parameter fitting and these coef- 
ficients can be exactly derived analytically. To summa- 
rize, there are two di fferences in the analy sis procedure 
from the Vela paper l|Enomoto et aLlEoOf}) : 

• events with hits in the outermost layer were re- 
jected ("edge cut"), and 

• OFF-source events were used for emulating the 
background behavior of the FD. 

In contrast, in the previous analysis, the background 
function was made from the ON data (outside of the 
source region), because of the lack of statistics of the 
OFF data. The FD has a position dependence and it 
deforms near the edge of camera. Thus with the high 
statistics OFF events, the current background determi- 
nation is improved and its position dependence inside the 
FOV of the camera was taken into account. The result of 
the fitting is shown by the blue points (the background- 
subtracted signal) and the red histogram (the best-fit sig- 
nal) in Fig.^ while the green is the best-fit background. 
We have a statistically significant signal within one de- 
gree of the SNR center. Note that this fitting procedure 
gives a maximum acceptance with reasonably small sta- 
tistical errors without introducing any image cuts, i.e., a 
fit gives the maximum acceptance of events remaining 
after the pre-selections and automatically minimizes the 
subtraction ranges. In addition, it removes any subjec- 
tive bias in cut determination. 

This method was checked by analysis of Crab nebula 
data taken in November 2004. The wobble-mode obser- 
vation was also used. The analyzable data corresponded 
to 316.4 min. The 9"^ plot is shown in Fig. 13 The number 
of excess events is 119 ± 22 (5.4cr). The flux is 1.2±0.3 
times the standard Crab flux with the power-law consis- 
tent with the standard index of —2.5. 

H.E.S.S. produced a very clean gamma-ray count map 
showing the shell-like structure of RX J0852.0-4622. In 



Fig. 4. — Morphology of gamma-ray-like events. The smoothing 
was carried out using the average of the center and neighboring 
eight pixel where pixel size was 0.2 X 0.2 degree^. The vertical scale 
(number of events per pixel) is indicated in the top bar. The dotted 
circle of 1 degree from the supernova center is shown. The dashed 
circle of 0.23 degree (plotted in lower left) is the (Icr) point spread 
function , where 68% of events are contained. The cross indicates 
the average pointing position, i.e., the center of the remnant and 
the squares the "wobble" pointing positions. The contours are the 
20, 45, 65, a nd 80% levels in the ASCA Gas Imaging Spectrometer 
X-ray map ITsunemi et aLEOOft : iSlane et alji2001iV 



making a 2D-count map, we tried a simple time-based 
subtraction. The normalization factor for OFF-source 
run is 1081/1031=1.05. We made two 2D-plots for arrival 
directions for ON and OFF using selected events with the 
FD in the range between —1 and 1, which is the signal 
region as shown in Fig. ^ The difference between two 
2D distributions with the above normalization is shown 
in Fig. 01 Comparing it with Fig. 2 of Aharonian et a^ 
(200.5^1), we confirm the shell- like structure, with less 
emission in the center position. Our data show that the 
emission in the northwest-rim (NW-rim) is the strongest. 
The signal-to-noise ratio (S/N) of our analysis is signif- 
icantly worse than that of H. E.S.S., which can be seen 
comparing the S/N of Fig. 1 inlAh aronian "eTaLl Il2005bfl 
and Fig. ^of this paper. Therefore, the statistical error 
in this map is larger than that in the H.E.S.S. map. The 
difference between red and green of this color map (Fig. 
0J) is not significant. 

Fig. 0] was made by positional subtraction. Fig. ^ was 
made based on fitting the event shapes (FD). We can 
also study the morphology by extracting gamma-ray-like 
events using the fitting procedure described so far. The 
FD distributions for each direction are recorded in mul- 
tiple histograms and the fittings were carried out. The 
resulting morphology is shown in Fig. |31 Although we 
observe a shell-like structure, the brightest region has 
changed from that of Fig. 21 As has been described, the 
difference between this figure and Fig. 21 can be consid- 
ered to be within the level of statistical fluctuations. 

We have thus confirmed the H.E.S.S. results that RX 
J0852. 0—4622 is a strong, shell-like TeV gamma-ray 
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Fig. 6. — Wide-range 6^ plot. The vertical scales are the num- 
bers of excess events per 0.3 degree^. The dashed histogram is the 
acceptance (in arbitrary units). 

Fig. 5. — Morphology of gamma-ray— like events obtained by the 
fitting procedure described in the text. The other definitions are 
as same as those of Fig. |1] 



source. We have also confirmed our previous claim of 
the CANGAROO-II detection of gamma-rays from the 
bright NW-rim (though not the reported soft spectrum) . 
The threshold of this analysis (Figs. ^ and ^ is esti- 
mated to be 1.02 TeV. This high threshold can be ex- 
plained by the deterioration of "light collection efficien- 
cies". The relatively poor S/N, even for the three- fold 
coincidence data, can be explained by the blur spot size 
of >0.1 degree of the segmented mirrors, which are made 
of plastic. Monte-Carlo simulations with zero spot size 
show a much better separation between gamma-ray and 
proton events. This means that replacing these mir- 
rors with glass or metal, we have a chance to approach 
H.E.S.S. quality even with 10 m diameter structure of the 
present telescopes. 

5. RESULT 

The radial distribution of the gamma-ray signal with 
respect to the SNR center is shown in Fig. 1 of 
lAharonian et alJ l|2005b) . We took the FD distributions 
of annuli of various squared-radii with respect to the SNR 
center and fitted them. The 0^ distribution over a wide 
range is shown in Fig. |B| The spatial excess appears only 
within one degree from the center of the SNR, in agree- 
ment with the H.E.S.S. result. The number of excess 
events is 557 ± 77 (7.2ct). 

The gamma-ray fiuxes within a one-degree radius were 
derived by fitting the FD distributions within this area 
on an energy bin- by-bin basis. The black squares in 
Fig. are obtained by this experiment. The numeri- 
cal data are listed in Tabled The consideration on the 
experimental energy resolution was carried out unfolding 
the generated energy from the detected photon-electron 
yield using the Monte-Carlo events under the assump- 
tion that the energy spectrum is proportional to 
(H.E.S.S. 's spectrum). The "Mean Energy" quoted in 
Table ^ is the mean generated energy of the accepted 
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Fig. 7. — Differential energy spectra; the red points by 
H.E.S.S. are for the whole remnant and the black points from these 
CANGAROO-III observations are also for the whole remnant. The 
error bars are statistical. 



TABLE 1 

Differential flux within one degree from the SNR center. 



Energy binning (TeV) Mean Energy (TeV) dF/dE (cm s' 



0.83 - 1.07 
1.07 - 1.37 
1.37 - 1.62 
1.62 - 1.91 
1.91 - 2.66 
2.66 - 3.70 
3.70 - 6.07 



1.02 
1.24 
1.51 
1.77 
2.24 
2.98 
4.72 



-^TeV- 
) xlO" 



( 2.77 ± 0.88) xlO-i^ 

( 1.32 ± 0.43) X 10-11 

( 1.03 ± 0.29) xlO"" 

( 4.40 ± 2.14) X 10-12 

( 7.11 ± 1.17) xlO-12 

( 1.93 ± 0.46) xlO-12 

( 8.83 ill 5.03) xlO-" 
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5JJS5 53395 53400 53405 53410 53415 

MJ.D (day) 

Fig. 8. — Time variability of the integral fluxes at the en- 
ergy greater than 1.02 TeV. The vertical axis is the flux nor- 
malized to the H.E.S.S. mean flux (dF/dE = 2.1 X IQ-^E-^ l 
[cm~^s~^TeV~^]). The horizontal axis is modified Julian days 
(M.J.D.). 



events in the Monte-Carlo simulation in each binning. 
The obtained differential flux is calculated at this " Mean 
Energy" under the assumption of the above power-law 
spectrum. The best fit with a power-law assumption for 
this energy spectrum is: 

dF 

— = [2.5 ± 0.6{stat.) ± 0.6(sys.)] x 10"" 

d-bj 



—) 

1 TeVy 



2.2±0.3{stat.)±Q.3{sys.) 



^TeV" 



Presently we estimated the systematic errors as follows; 
the main one is the energy scale of 15% due to the abso- 
lute light collection efficiency and the rest is the detector 
systematics of 10%. The ambiguities of the absolute flux 
and the power- law strongly corre l ated to each other. The 
data points of lAharonian et all l)2005bD were obtained 
from 

http : //www. mpi — hd. mpg.de/hfm/HESS/HESS. html 

Our results are consistent with H.E.S.S. for both flux 
level and spectrum, i.e., the CANGAROO-III results for 
the whole remnant yield a harder spectru r n tha n the 
CANGAROO-H spectrum of IKatagiri erall 

We checked the time variability of these excesses. Five 
periods were selected so that they have similar statistics, 
i.e., two-to-four days periods. The integral fluxes at en- 
ergy greater than 1.02 TeV, which are normalized to the 
H.E.S.S. mean flux are plotted in Fig. |S| It is consistent 
with stable emission. 

We now consider the fluxes and spectra as a function 
of position in the SNR to attempt to reconcile the dif- 
ferences between the CANGAROO-H and H.E.S.S. re- 
sults. The energy spectrum obtained by CANGAROO-H 
is shown by the red points in Fig. The CANGAROO- 
H flux is barely consistent with the other results. The 
only difference is that the single telescope CANGAROO- 
H observations yielded an excess only at the NW-rim. 
The exact procedure is reviewed as follows; 

• Derive Hillas parameters and calculate 
Likelihood Ratio (LR). Probability Density 
Functions (PDFs) were derived for both gamma- 
ray and cosmic ray initiated events. The PDFs 
for gamma-rays were obtained from simulations, 
while those for cosmic rays were obtained from 
OFF-source data, both in energy-by-energy 
bases. Histograms were made of length and 
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Fig. 9. — Differential energy spectra; the red points are obtained 
by CANGAROO-II around the NW-rim and the black points are 
obtained by this observation for the quadrant centered around the 
NW-rim. The yellow area is 1/4 of the H.E.S.S. fiux. 



width using both data sets; these distributions 
were then normalized to unity. The probability 
(L) for each assumption was thus obtained by 
muhiplying FBF{width) by PBFilength). In 
order to obtain a single parameter, and also to 
normalize it to unity, we used the Likelihood- 
ratio (LR): -r? — ; r / t — r, therefore, m 

V / L{gamma—ray)-\-L{proton) ' ' 

energy-by-energy bases. 

• Select LR > 0.4 events and plot Alpha{"lma.ge ori- 
ented angle"). 

• Normalize ON- and OFF-source distributions at 
Alpha > 30 degrees and obtain the excess within 
Alpha < 30 degree. We did not use the time- 
based normalization (considering the ON- and 
OFF-source observation time ratios). Only the 
shape of the OFF-source Alpha distribution was 
used to determine the normalization factor. 

Only in the case that the spectral index in the NW-rim is 
softer than elsewhere can the softer spectrum obtained 
by CANGAROO-II be explained. We counted the ex- 
cess within the quadrant within ±45 degree from the 
NW-rim center (RA, dec = 132.245°, -45.650° J2000). 
The result is shown by the black points in Fig. O which 
is slightly higher than 1/4 of the whole emission (the 
hatched area in Fig. O . Its spectrum is still consistent 
with that of the whole remnant. Therefore, it is not 
consistent with the CANGAROO-II flux within the sta- 
tistical errors. Additional sys t emati c uncertainties such 
as outlined in IKatagiri et alJ l)2005f) are required. The 
main one was considered to be energy uncertainty due to 
the deterioration of mirrors and PMTs. Also saturation 
of the electronics might contribute. In particular, the 
lower energy fluxes were derived with large acceptance 
correction factors and have larger systematic errors. 

We analyzed the same data again using a likelihood- 
ratio instead of the Fisher Discriminant for the 
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Fig. 10. — Wide range 9^ plot obtained using the "Likelihood 
analysis" . 
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Fig. 11. — Differential energy spectra; the red points were ob- 
tained using the "Likelihood- Ratio analysis" for the whole remnant 
and the black points by the FD analysis for the the whole remnant. 
The error bars are statistical. 



gamma/hadron separation in order to check the 
CANGAROO-II analysis. Here we used lengths and 
widths of three telescopes and used a safety cut of 
LR > 0.1 for this supernova remnant in Fig. ^1 We 
see a similar excess within the remnant. The number of 
excess events is 538 ± 128 (4.2(7), consistent with that 
obtained by the FD analysis. The flux obtained by this 
method is shown in Fig. The best fit with a power- 
law assumption is: ^ — [2.4 ± 0.8{stat.) ± 0.6{sys.)] x 



10- 



(l Tcv) 



2.0±0 A{stat.)±0.3{sys.) 



Therefore, the detection of this supernova remnant by 
iKatagiri ct al. (200^ is shown to be correct. Hard- 
ware correction for the "light collecting" efficiencies 
and "saturation" effect have been refined since the 
CANGAROO-H resuhs, and since the CANGAROO-H 



telescope (Tl) has not been operated in recent years, 
it is hard to estimate these deteriorations over 2001 to 
2003. 

6. DISCUSSION 

We now discuss the possible origins of these TeV 
gamma rays. The key constraints are the energy of the 
protonic cosmic rays are ~ 10^" and the electronic ~ 10^* 
erg. The former is required for the supernova origin of 
cosmic rays and the ratio between the former and lat- 
ter is b ased on the electron m easurements in the solar 
system (iKobavashi et alJl2004|) . We, here, argue the dis- 
tance of this SNR and the mechanism of the gamma-ray 
emission. 

We first modeled the TeV gamma-ray flux assum- 
ing Inverse Compton scattering on the cosmic mi- 
crowave background. The input electron spectrum is 
assumed to be proportional to E~'^ e~^^^'"''''' . Here 
the number of the data points are 16 (H.E.S.S.=9 and 
CANGAR00-HI=7). The best-fit results under a fixed 
parameter of 7 = 2.1 (a typical assumption) are for 
an exponential cut-off energy of Emax = 26 ± 5 TeV 
with x^/o?o/=13.9/(16— 2). Under the assumption of 
d=200pc (distance of this SNR from earth) we ob- 
tained total electron energy at E greater than 1 GeV 
of JE^e = 0.04 ± 0.01 • 10'*^ erg which is low considering 
the supernova origin of cosmic rays under the assump- 
tion of the elec tron/p roton ratio (e/p) at the acceleration 
site l|R,evnoldjll996t IKob avashi et alJl2004D . Note that 
we only consider the statistical errors of the observation. 
The real errors of the fitted parameters would be larger 
than those cited. 

Second we consider protonic contributions to the TeV 
flux. Assuming a shock velocity of 5000 km s^^ in the 
free expansion era and a swept mass of order Mq with 
the age of this SNR of 630 yr (d=200pc), the ambi- 
ent interstellar matter (ISM) density of n=:0.2p/cc is 
obtained. In the case of d=lkpc, n=0.002p/cc. We 
carried out a two-parameter fit {Emax and Ep) to only 
the TeV flux with a fixed parameter of 7 — 2.1 where 
the input proton flux was assumed again to be propor- 
tional to £;-Te~-^/-^'""== . The best fit was obtained for 
Emax > 70 TeV with x^/dof = 15.2/(16 - 2). The pro- 
tonic CR energy at d=200pc is Ep = 0.80± 0.08 • 10^^° erg 
under the assumption of the ambient ISM density of 
0.2 p/cc at d=200pc. We can not reject a protonic contri- 
bution to the TeV emission for d=200pc. We, however, 
should have a significant electronic contribution and to 
produce a significant protonic component we need to sup- 
press the e/p-ratio to far less than 0.01. When we con- 
strain the total cosmic-ray energy to be less than 10^" erg 
at d=lkpc, the ambient density of the ISM is required 



Ip/cc 



M, 



o 



to be at least 4 p/cc. A total mass of '--^ 550 
within a 1 degree radius is necessary in that volume with 
a shell-like distribution. This kind of heavy of molecu- 
lar cloud is, however, marginally within the detectable 
range of recent CO observations. The search for such a 
molecular cloud is necessary. 

We now shift to the lower energy region. The 
radio data used here are based on Table 2 of 
iDuncan and Green ( 20M) and the X-ray data are from 
Table 1 o f Slane et idl l|20?)l (ASCA). Note that 
iSlane et"al.i (2001) show fluxes for the north-west, north- 
east, and west regions of this SNR. We combined them 
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Fig. 12. — Spectral energy distributions. The data points are 
for radio, X-ray, H.E.S.S., and this observation. The colored curves 
are model calculations; the blue at low energy are for synchrotron 
radiation and at high energy are Inverse Compton, the red are -k^ 
decay, the green, Bremsstrahlung, and the light blue show the sum 
at high energy. 



and calculate five data points from 0.5 to lOkeV. They 
only gave errors for the spectral indices. We, therefore, 
take only them into account, i.e., the errors are under- 
estimated. With the above best-fit parameters under 
the electron assumption, an ambient magnetic field of 
B =~ 4 also gives an allowable one-parameter fit to 
the X-ray spectrum and lower estimation to the radio 
observation. 

A filament al X-ray structure w as observed for this SNR 
by Chandra fBam ba et all l2005y>. We might need to in- 
troduce a magnetic field filling factor (/b < 1). The 
best-fit results freeing all five parameters are shown in 
Fig. El The ^jdoj was 18.8/(23-5) with 23 data 
points. The fitted spectral energy distribution is shown 
in Fig. El The blue lines are the electronic contribu- 
tions described above. We obtained: B = 5.8 ± 0.3 /iG, 
!b = 0.40 ± 0.04, 7 = 2.37 ± 0.01, E^^^ = 37 ± 2 TeV, 
and = 0.18 ± 0.02 x IC** erg at d=200pc, less ener- 
getic than that of the standard supernova origin. Also 
the filling factor looks too large to create the X-ray fila- 
mental structure. If this SNR is assumed to be located 
at d=lkpc, the energy simply becomes 25 times bigger, 
that is, marginally higher than usually assumed. This 
low magnetic field strength with the estimated age of 
630 yr at d=200pc is considered to be within the al- 
lowable range for theory of electron acceleration up to 
40 TeV. Assuming the proton spectrum is the same as 
that of the electrons and its energy is hundred times 
larger, the red line is the expectation of the 7r° — > 77 con- 
tribution which is a negligibly small contribution to the 
TeV emission. The green line is that for Bremsstrahlung. 
Note that this calculation has a big uncertainty in the 
sub-MeV region. Suzaku observations of the hard-X re- 
gion are eagerly awaited. In the case of 7 = 2.1, the 
contribution from Bremsstrahlung is below the Suzaku 
sensitivity limit. 

Conside ring the filamental structure observed in X- 
rays, Bamba et al.l l)2005t) assumed a magnetic field as 
high as 500 /xG. This model predicted the maximum en- 
ergy of the electrons to be Emax = 3 TeV. With that 
energy spectrum, the observed flat TeV emission can not 
be well described via either Inverse Compton or vr*^ de- 
cay processes. In order to fit the TeV spectrum, an or- 
der of magnitude higher cut-off energy is definitely nec- 



essary, which favors a two orders of magnitude lower 
magnetic field, because the maximum energy of the syn- 
chrotron radiation is a function of EmaxV^- This con- 
tradiction can be resolved by introducing a different en- 
ergy spectrum between electrons and protons as was pre- 
dicted by kinetic non-linear acceleration theories where 
synchrotron cooling of electrons was tak en into account 
(see for example iBerezhko et al. | l|2003j) '). The filamen- 
tal structure might be a result of strong synchrotron 
losses of emitti ng electrons in am plified fields as i n case 
of other SNRs llVink fc Laming 1.2003: Berezhko et al. I 
120031 IBerezhko fc Volk 1 120041) Such a high magnetic 
field (hundreds of /iG) in the downstream region reject 
an Inverse Compton solution thus supporting the idea of 
a nucleonic origin of the observed TeV gamma-ray flux. 
However, in that case, we need an origin other than SNRs 
for the observed electrons in the solar s ystem. Our sim- 
ple es timation agrees with their models l|Berezhko et al. I 

Although at present a perfect interpretation is missing, 
fine angular resolution studies over a broad energy range 
should help to determine one. 

7. CONCLUSION 

We have observed the supernova remnant RX 
J0852.0-4622 with the CANGAROO-III stereoscopic 
imaging Cherenkov telescope. It is clear from the results 
for RX J0852.0-4622 that CANGAROO-III's stereo- 
scopic performance is significantly improved from that 
of the single CANGAROO-II telescope. We have con- 
firmed the H.E.S.S. group's result about the morphology 
and the energy spectrum. The NW rim is the bright- 
est feature. RX J0852. 0—4622 has been confirmed to 
be one of the strongest TeV gamma-ray sources in the 
southern hemisphere and shows a clear shell-type struc- 
ture at these energies. The energy spectrum around the 
NW-rim was measured to be consistent with that of the 
whole remnant, not in good agreement with the previous 
CANGAROO-II result. The difference can be partially 
explained by the deterioration of the hardware of the 
CANGAROO-II telescope. 
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APPENDIX: 2004 OBSERVATION RESULTS 

Additional observations were carried out in the period 
from 2004 January 17 to February 25 in the same "wobble 
mode" as the 2005 observations. The average pointing 
position was the NW-rim (RA, dec = 132.245°, -45.650° 
J2000), i.e., the maximum of the X-ray emission. How- 
ever, only the second and third telescopes, T2 and T3, 
were used for the 2004 observations as T4 became opera- 
tional after these observations had been completed. Also 
in this period the "Global Trigger" system was not ready. 
Data were recorded for T2 or T3 independently when 
the photomultiplier signals exceeded 7.6 photoelectrons 
(p.e.) in more than four camera pixels ( "Local Trigger" ) . 
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The GPS time stamp was recorded in each dataset. An 
offline co incidence of time starn ps from T2 and T3 within 
±100 /iS l|Enomoto et alJl2006D was required for a stereo 
event. The typical trigger rate for each telescope was 
80 Hz, which was reduced to 12 Hz for stereo events. The 
light-collecting efficiency was calibrated during the offline 
processing by analyzing /x-rings and found to be 70±5% 
for both T2 and T3, a s for the Vela Pulsar observations 
ijEnomoto et aljf2006j) . The total observation time was 
1781, and 436 min, for ON- and OFF-source data, respec- 
tively. After pre-selection, the event rate was reduced to 
10 Hz. Finally after the "cloud and elevation cut" the 
analyzable data (ON-source run) corresponded to 1120 
min. with a mean elevation angle of 70.3°. The OFF- 
source data, however, only had 379 min., in part as these 
observations were carried out at a similar time with the 
CANGARO O-HI observations of t he Vela pulsar region 
described by lEnomoto et all l)2006j) . 

The differences between these and 2005 data are sum- 
marized as follows; 

• 2-fold (software) coincidence in 2004 and 3-fold 
(hardware) coincidence in 2005, and 

• the OFF-source data in 2004 are only 1/3 of the 
ON-source data. 

Although the light collecting efficiency is higher than 
2005, there remain some uncertainties, especially in the 
background subtraction. We, therefore, report these re- 
sults in this Appendix as supporting evidence to the main 



arguments. 

The same kind of analysis as for the 2005 data were 
carried out. Here, we restrict the analysis only to the 
higher acceptance regions due to the lack of statistics 
in OFF-source runs. The higher acceptance region was 
chosen as within 1.4 degrees from each "wobble" point- 
ing. Only the overlapping region was analyzed. This is 
approximately the north-east half of the remnant. Com- 
paring with the 3-fold data, the signal-to-noise ratio is 
poor. The 3-fold coincidence data had better background 
rejection power than that of the 2-fold data. The errors 
in the subtracted data are dominated by the OFF-source 
statistical errors. Due to these, the statistically accurate 
background functions could not be obtained position-by- 
position bases in the FOV. The background also behaved 
non-uniformly near the edge of the camera. Therefore, 
we did not use the fitting procedure of the 2005 anal- 
ysis. The ON— OFF subtraction was carried out using 
time-based normalization factor. 

As has been described, we can only analyze half of the 
remnant. The flux inside the maximum acceptance re- 
gion is shown in Fig. 1131 i.e., fluxes of the half remnant. 
The energy resolution was estimated using Monte-Carlo 
simulations in the same way as the 2005 analysis to be 
AS/i; = 17/^£;/l TeV+12%, which is slightly worse 
than the 3-fold analysis. The energy spectrum is simi- 
lar and the flux is half compared to the whole remnant. 
This is supporting evidence for the results of the 2005 
observations. 
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Fig. 13. — Differential energy spectra; the black points from the 
3-fold coincidence analysis of 2005 CANGAROO-III data for the 
whole remnant (the same as the black points in Fig.|2l and the red 
points from the 2-fold analysis of 2004 CANGAROO-III for about 
a half of the remnant. 



